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Abstract 



The low-energy properties of the Anderson impurity are studied under a finite 
bias voltage V using the perturbation theory in U of Yamada and Yosida 
in the nonequilibrium Keldysh diagrammatic formalism. The self-energy is 
calculated exactly up to terms of order uj 2 , T 2 and V 2 using Ward identities. 
The coefficients are defined with respect to the equilibrium ground state, and 
contain all contributions of the perturbation series. From these results, the 
nonlinear response of the current through the impurity has been deduced up 
to order V s . 
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The Kondo effect |T| in quantum dots is a very active field of current research, and 
the results of recent experiments 0] have been shown to be in qualitative agreement with 
theoretical predictions ||. There is a new feature, the Kondo effect under a finite bias voltage 
V, which is specific to quantum dots and has no analogue in dilute magnetic alloys. So far, 
a number of theoretical methods have been used to investigate out of equilibrium Kondo 
physics: perturbation theory in the Coulomb interaction U of the Anderson model 0], the 
noncrossing approximation in the limit of U — > oo 0], the bosonization approach based on 
the Kondo model 0, and so on. Although these works give some important insights into 
the electron correlation in a driven system, further information is needed to clarify these 
new aspects of Kondo physics. 

In the equilibrium state V = and linear response for V ^ 0, the low temperature 
properties of the Kondo physics can be described by a local Fermi liquid theory [[7|,|8|. The 
microscopic basis for this theory has been provided by a diagrammatic analysis based on 
the assumption that the ground state evolves continuously with the adiabatic switching-on 
of the Coulomb interaction U ||||. Some of the features found from the first few terms of 
the perturbation series have been shown to hold to all orders in U, and can be expressed 
in terms of Fermi-liquid parameters which characterize the quasi-particle excitations. The 
perturbation theory results have been confirmed by the exact Bethe ansatz solution jlU 



The local Fermi liquid theory also gives a description of the linear response of the current 



through small interacting systems fill]. In the nonlinear response, Hershfield et al have given 
the asymptotic form of the self-energy of the Anderson model to the second-order in U for 
the electron- hole symmetric case 0. In this paper, we show that some of the features of 
the nonequilibrium state hold to all orders in U without the assumption of the electron-hole 
symmetry. Our proof uses Ward identities for the derivative of the self-energy with respect to 
V, derived within a diagrammatic analysis based on the Keldysh formalism. Consequently, 
the low-energy behavior of the self-energy £ r (cu) can be expressed in terms of Fermi-liquid 
parameters for the equilibrium ground state; Eqs. (|Hf)-(|T§|). The nonlinear response of the 
current J has been calculated up to order V 3 . In the electron-hole symmetric case, the 
expressions for £ r (u;) and J can be written in simplified forms Eqs. (|20| ) - (|2ll) . 

We start with the Anderson impurity connected to two reservoirs at the left (L) and 
right (R); 

H = - J2 4 c \<r c ja + E J2n a + U n^n ol 

\=L,R ijeX cr 
a 

~J2 V L [4 C -l<x + cl la C 0a ] - J2 V R [4 C 0<7 + 4^ } • (1) 

Here tfj (t^j) is the hopping matrix element in the left (right) reservoir, c| CT is creation 
operator of an electron with spin a at site i, and n 0o . = CQ a c 0cr . The couplings between the 
impurity and reservoirs are described by vl and vr. The site indices are assigned to be 
i < —1 for i G L, and i > 1 for i G R. The potential due to the bias voltage V is included 
in the diagonal element t^. We assume that Eq is a constant independent of V, and take 
the chemical potential in the equilibrium state, /x, to be the origin of the energy. 

The Keldysh formalism is described by the four types of Green's functions at the 
impurity site 0: G~ (t) = —i (T c 0a (t) Cq ct (0)), G~+(t) = i (4(0) c 0(7 (t)>, G+~(t) = 
—i (c 0a (t) 4(0)), G ++ (t) = —i (Tc 0cr (t) 4(0)), and c 0a (t) is the Heisenberg operator which 
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at t = coincides with the interaction representation. Furthermore, the retarded and ad- 
vanced functions are given by G r = G — G~ + and G a = G — G H . The average (• • •) is 
taken with the density matrix at t = 0. Initially at t = — oo the reservoirs are separated from 
the impurity and described by the thermal equilibrium of the chemical potential //£, and /i R 
with fii — Hn = eV. Then vl, v r and U are switched on adiabatically. In the noninteracting 
case U = 0, the Fourier transform of the retarded and lesser Green's functions are given by 



G 



G r {u J ) = [u J -E + i{T L + T R )} 
f L {uj)T L + f R {uj)T R 



-i 



LU = - 



r L + r 



n 



R 



G a (u) 



(2) 
(3) 



Here /^(w) = f(u - fi LfR ) and f(u) = [e w/ +1] . T L = np L v L , T R = np R v R , and 
Pl (Pr) 15 the local density of states at the interface i — — 1 (+1). We assume p LR is a 
constant and the bandwidth is infinity, so that is a constant independent of u. From 
Eqs. (0)-(0), other noninteracting Green's functions can be obtained using the properties 
G~ + + G+- = G + C++, G ++ (lu) = -{G—(u)}*, and G a (u) = {G r (u)}*. Note that the 
distribution function is introduced through the fiXi + fR^R dependence of Eq. (|]). It is not 
a symmetric function of pi and p R if Y L ^ T R , and thus the local charge at the impurity 
site is affected by the inversion of the bias voltage. The full Green's functions are described 
by the Dyson equation [Giu)}^ 1 = {Gq(co , )} _1 — S(u;) ; 



Go 



G G + 
c<-\ — 

Lr Ltq 



E++ 



(4) 



Here S(u;) is the self-energy due to U, which can be defined using the diagrammatic repre- 
sentation of the perturbation series P|,l2|. Furthermore, the retarded function is given by 
{G r } _1 = {Go} -1 — E r with E r = E + E~ + , and the current flowing through the impurity 
can be written in terms of the spectral function A(u) = —(1/tt) ImG r (w) [13]; 



R 



J= 2e 4T L : 
h T R + T L 



du[f L {Lu) - f R (u)]irA(u). 



(5) 



Generally A(u) is not symmetric against the inversion of the bias voltage due to the fiXL + 
f R T R dependence of Gq v , and thus the nonlinear part of J transforms to —J only for 
Y L = T R . In order to specify the deviations of pl and p R from the equilibrium value, we 
introduce the parameters qj, and a R : p L = ai eV and p R = —a R eV with a>L + a R = 1. The 
inversion of the bias voltage is described by the transformation (V, q^, a R ) =3- (—V, a R , o>l). 
In the following, we will discuss the behavior of E r (u;) in the vicinity of the equilibrium up 
to order V 2 by evaluating the first and second derivatives with respect to V. 
We now consider the first derivative. At V = 0, Eq. 



Gn 



v=o 



-m 



G r 0eq (u) 



G a 0eg (u;) 



With Goegi^) 



= G5( W ) 



is written as G ^ q {uo) = 
The derivative of Gq u ' 



with respect to V can be written in terms of the Green's functions for the equilibrium state; 




v=o 




dE r oeq 



Goeqi^) 



Go eq {u) 



(6) 
(7) 
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Here a = (a L T L -a R T R )/(T L + T R ) and G t 



Oeg 



UJ 



Cr 



{ ■X ) . with v, v' = +, — . Note that 
The derivative of XI (a;) can 
taking the derivative of Gq 



v=o 

T L /(T L + T R ). 

q u 's within 

In the closed loops of the diagrams, the frequency of Gq u (uii) 



a becomes zero for oil = T r /(Tl + T r ) and a R 
be carried out using the diagrammatic analysis, i.e 
the diagrams for S(a>) 0. In the closed loops of t±±c u. 
can be shifted to uii + li without changing the value of E(a>) because the frequencies of the 
closed loops are integrated out. Therefore, it is possible to assign the frequencies of the 
diagrams so that all the noninteracting Green's functions possess the external frequency ui 
in the argument. Then, the derivative d/d(eV) can be replaced by —a (d/dui + d/dE ) for 
all the internal Gq u (uii + w)'s using Eq. (§), and thus 



<9£M 



d(eV) 



— a 



v=o 




T, eq (ui). 



(8) 



Here Y, eq {uj) 



= Si 



\v=o- 



We will assign the label 'eg' to the subscript of the functions for 



the equilibrium state. Note that Eq. (§) holds at any temperatures since no assumption has 
been made for the value of T so far. The same relation holds for E r (cc;), and operating d/du 
to the expression corresponding to Eq. we obtain 



<9 2 £ r 



.oj 



dcud(eV) 



v=o 




In the equilibrium state, the imaginary part of the self-energy behaves as ImS 



b 2 + 



(ttT) 2 ] for small u and T §g|TJ 



eq 



(9) 



Thus, substituting this into the right-hand-side of 
Eq. (^|), we find that ImE r (w) has the cross-term ueV. Similarly, differentiating Eq. (§) 
with respect to T, we find that there is no cross-term of TeV in the imaginary part. At 
T = and V = 0, the usual diagram technique in terms of G~~ is applicable, and the causal 
element of Eq. (|8|) can be written in terms of the vertex part using Eq. (0) |§, 

( d 




eq 



( UJ 



E 



duo 1 
2iri 



Faa';a>a{uJ, Ul' \ Ul' ', Ul) 8G{u') 



df(ui') 
dui' 



5G(0) 
2-ni 



a(T . .yo-iui, 0;0, uj) 



(10) 



Here r cr(T / ;o -/ cr (ct;, ui'\ uj', uj) is the antisymmetrized total vertex part of the T = formalism, 
5G(uj) = G' eg (ui) - G a eq (ui), and -df(u)/du = 5(uj). We note that at T = the differential 
coefficients are related to the enhancement factors of the physical quantities: 7 = 1 — 
ui)/duj\ u;= o 1 Xe = XTT + Xn = 1 + d K q (Q)/ dE o, and Xs = XtT ~ Xu> where X™' = 

with h a i being an external field described by H ex = —Y> a h a n§ a . 



<W - dEl q>a (0)/dK 



Furthermore, there are some exact relations among these parameters [§,|§: Xtt = 7 an( ^ 
X U = -A eg (0)r U;iT (0,0;0,0). 

Next we consider the second derivative with respect to V in the similar way. The second 
derivative of the noninteracting Green's functions can be written as 




K 



v=o 



' d_ _d_ 

dui dE 



Goeqi^) 



GZ^ui) 



d 2 f(uj) 
dui 2 



G r 0eq (uj) 



G a 0eq (uj) 



(11) 
(12) 
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Here k = {ol 2 l Yl + ^r^rJ/O^l + Fr). The second derivative of S(u;) can be evaluated by 
operating {d/d(eV)} 2 to Gq u ' {u>i+u>ys within the diagrams for £(u;). The contributions are 
classified into two groups according to whether {d/d(eV)} 2 operates to, (i) two different 
electron lines, or (ii) a single electron line. In the class (i), each of the first derivatives 
d/d(eV) can be replaced by —a{d/du + d/dE ). In the class (ii), the second derivative 
d 2 /d(eV) 2 can be replaced by k (d/dtu + d/dE ) 2 . Therefore, 



d 2 Y,{uj) 



d(eVy 



v=o 



-0, 

a'id + Cn) + (K-a 2 )C n 

2 



a 



' d_ 







du dE, 







(r L + r R 



2 D 2 ^^) 



(13) 



Here a 2 Cj and k C n denote the contributions of the classes (i) and (ii), respectively. Note 



that Cj + C n = (d/du + d/dE ) 



S e q(cu) 



Cu = D 2 ~S eq (uj), and D 2 denotes the functional 



<■<! 



operation carrying out (d/du + d/dE ) for all the single Gq u . At T 
zero-temperature formalism is applicable for the causal element Z} 2 £ 
operation to pick up one electron line from the diagrams for E 
Eq. p|. 
obtain 



eq 



Carrying out the same operation and using Eq. (|T 



XI / ^^aa';a'a(uJ,Uj';Uj',Uj) 5G(uj') I 
a' K 



= 0, the usual 
The functional 
has been used to derive 
instead of Eq. (0), we 



5G(0) 9 , 




w'=0 



2ni 



^O A ( W) 0;0,u;). (14) 



Here SG'(u) 



I" ' a ' a\aa' , UJ', UJ , Cj' 



d5G(u))/duj, and —d 2 f{uj)/d(jj 2 = S'(u). Note that r CTCr / ;(T / CT (c<;, cj'; cj', uj) = 
due to the antisymmetric property of the total vertex. Furthermore, 
owing to the rotational symmetry of the spins, the first term of Eq. ([14]) can also be written as 
J2cr'{9/duj') T aa i. a i a (uj f ,uj;uj,uj'). Now we consider the limit to — > 0. As shown by Eliashberg 
quit generally, the vertex part T aa '. a / a (uj, uj'; uj', uj) has singularities at ui—ui' = and ui+ui' = 
||15|| . In Fig. [15], the diagrams contribute to the singularities for small u and uj' are 
shown: (a)-(b) and (c) yield the \ui — ui'\ and \ui + ui'\ dependence, respectively. In the 
figures, the solid line denotes the full Green's function, and the square denotes the vertex 
part with zero frequencies ^.^(O, 0; 0, 0). Note that r^.^O, 0; 0, 0) = because of the 
Pauli principle. Therefore, the derivative of the vertex part has a discontinuity of the form 
(d/du') \ui — uj'\ = — sgn(u — uj'), and the limit uj, uj' — ► depends on which frequency is 
first taken to be zero. These discontinuities appear in the imaginary part, and for small uj 
and uj' the contributions of the diagrams (a)-(c) are written as 

d 

X — Im r ua ,.y a (uj, uj'; uj', uj) 



x Im 



r U;iT (o,o ; 0,0)| 

du" 







G--{uj")—G--{uj-uj' + uj") + 



2ni ~ eq y " ' du' 
~ n Xu [ _2 s g n ( w ' - uj) + sgn(u/ + uj) 



fit i" r) 



(15) 
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Note that G~{u) 
first, we obtain 



ReG r eq (uj) +ilmG r eq (uj) sgn(uj) at T 



0. Thus, taking the limit uj' — > 



lim ImD 2 S m 



[u) = -3 7rA eg (0) x 2 u sgn(c^). 



(16) 



The opposite limit of Eq. ([15]) corresponds to the uj 2 term of ImE (uj), which is obtained 
by operating d/duj to Eq. (0) as (d 2 /duj 2 ) Im£--(a;)| = -7rA eg (0) Xn sgn(cu) g|. On 
the other hand, the real part of the vertex part does not show these discontinuities, and the 
limit is independent of the ways for uj, uj' — > 0. Thus, in the real part of Eq. (|I~4D , we can 
take the limit uj — > first. Then using Eq. (|10D , we obtain 



Re^E" 



d_(d_ _d_ 

U " duj [du + dE 



Re E 



eq 



[UJ 



uj=0 



A eg (0) 



xn- 



(17) 



Consequently, using Eqs. (|i3|) and (|I6|)-(|T7|), <9 2 S e(J /d(eV) 2 can be expressed in terms 
of the enhancement factors for the equilibrium ground state. 

We now summarize the results to show the low-energy behavior of E r (u;). Since the 
imaginary part has no TeV dependence as mentioned, the result, which is valid up to uj 2 , 
V 2 , and T 2 , can be written as 



nA eq (0) x 2 n 



' uj — a eV ) + 



ST L T R 



(r L + r 



(elO 2 + (nTf 



+ 



(18) 



Note that electron-hole symmetry has not been assumed so far. The general expression of 
the real part is rather long. At T = it takes the form 



Re E r M 



s:,(o) + 



eq\ 



~ \ ~ i, u 2 

X TT )^ - OiXneV + -u - a 



dxn 



ueV 



1 

+ 2 



a 



dXs 
dE 



Tr,T 



L-L R 



0E + A eq (0) XU 



[eV) 2 + ---. (19) 



Here b = (d 2 /duJ 2 )ReT,l q (uj) . In the electron-hole symmetric case (E = —U/2, Tl 
Tr, and = a R = 1/2), the real part Eq. ( |T9"D simplifies and we obtain the expression, 



{UJ 



uj 2 + °- (eV) 2 + (vrT) 2 



+ 



(20) 



which for V = corresponds to the result of Yamada and Yosida II]. Here A = Tr + Tl. We 
note that in the symmetric case the Bethe ansatz solution for x acr > can be expanded in power 



series of U which converges absolutely for any finite U ||10|| , and the second-order result of 
Im E r (w) of Hershfield et al is reproduced through the first term of the perturbation 
series xn = —U/{n/S) + • ■ ■ PJIT|. In the general formula Eqs. (|T8|)-(pJ|), the contributions 
of the higher order terms are described through the parameters x a a', A eq (0), YT e (0), and b. 

Using the results of the self-energy Eqs. ([T8"|)-(|T9"|), the low energy behavior of A{uj) is 
obtained correctly up to terms of order uj 2 and V 2 . Then, using Eq. (^j, the current is 
evaluated exactly up to order V 3 as J = g{V + g?V 2 + g^V 3 + ■ ■ -. Note that the current 
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conservation is fulfilled up to order V 3 because all contributions of the perturbation series 
are taken into account for the coefficients g%, #2 and g^. Generally, the explicit form of g% 
is rather long reflecting the lengthy expression of the real part Eq. (|19|). Specifically, in the 
electron-hole symmetric case, gi vanishes, and g\ and g% can be expressed in terms of the 
two parameters xii an d Xn, an d the resulting expression of the current is 



J 



h 



v 



i - 



X^ + 2 X 2 n (it_T_ 
A 



12 



A J 



+ 



(21) 



Furthermore, in the Kondo limit (U — > oo with Eq = —U/2), the charge fluctuation is 
suppressed \ c — > and the low-energy behavior of the current is characterized by the single 
parameter \ s as 



2e 2 



4 V A 



Xs f eV \ 



(22) 



In this limit the spin susceptibility is determined by the Kondo temperature as Xs A/T^, 
and thus the V 2 and T 2 dependence of the differential conductance dJ/dV is scaled by Tk- 

In conclusion, we have derived the Ward identities for the derivative of the self-energy 
with respect to V based on the perturbation theory in U. Using the identities, we have 
calculated exactly the low-energy behavior of the self-energy and nonlinear response of the 
current in the nonequilibrium steady state. 

I would like to thank A. C. Hewson for useful comments and discussions, and H. Ishii and 
S. Nonoyama for valuable discussions. I wish to thank the Newton Institute in Cambridge 
for hospitality during my stay on the programme 'Strongly Correlated Electron Systems'. 
This work is supported by the Grant-in-Aid for Scientific Research from the Ministry of 
Education, Science and Culture, Japan. 



7 



REFERENCES 



[1] A. C. Hewson, The Kondo Problem to Heavy Fermions (Cambridge University Press, 
Cambridge, 1993). 

[2] D. Goldharber-Gordon, H. Shtrikman, D. Mahalu, D. Abusch-Magder, U. Meirav, and 
M. A. Kastner, Nature 391, 156 (1998); S. M. Cronenwett, T. H. Oosterkamp, and L. 
P. Kouwenhoven, Science 281, 540 (1998); F. Simmel et al, Phys. Rev. Lett. 83, 804 
(1999). 

[3] T. K. Ng and P. A. Lee, Phys. Rev. Lett. 61, 1768 (1988); L. I. Glazman and M. E. 

Raikh, Pis'ma Zh. Eksp. Teor. Fiz. 47, 378 (1988) [JETP Lett. 47, 452 (1988)]. 
[4] S. Hersfield, J. H. Davies, and J. W. Wilkins, Phys. Rev. B 46, 7046 (1992). 
[5] N. S. Wingreen and Y. Meir, Phys. Rev. B 49, 11 040 (1994). 
[6] A. Schiller and S. Hersfield, Phys. Rev. B 51, 12 896 (1995). 
[7] P. Nozieres, J. Low Temp. Phys. 17, 31 (1974). 

[8] K. Yamada, Prog. Theor. Phys. 53, 970 (1975); Prog. Theor. Phys. 54, 316 (1975); K. 

Yosida and K. Yamada, Prog. Theor. Phys. 53, 1286 (1975). 
[9] A. Yoshimori, Prog. Theor. Phys. 55, 67 (1976). 
[10] V. Zlatic and V. Horvatic, Phys. Rev. B 28, 6904 (1983). 
[11] A. Oguri, Phys. Rev. B 63, 115305 (2001); J. Phys. Soc. Jpn. 66, 1427 (1997). 
[12] L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964) [Sov. Phys. JETP 20, 1018 (1965)]; 

C. Caroli, R. Combescot, P. Nozieres, and D. Saint- James: J. Phys. C 4, 916 (1971). We 

have used the notation given by E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics 

(Pergamon Press, Oxford, 1981). 
[13] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). Note that their definition 

of and Tpi are different slight from ours (by a factor 2). 
[14] The low-energy behavior of ImE^(w) given in the electron-hole symmetric case in || 

holds also in the asymmetric case just by replacing the coefficient xldd/ ^ by a general 

expression X||7rA e? (0) 0. See also discussions below Eq. flT5]). 
[15] G. M. Eliashberg, Zh. Eksp. Teor. Fiz. 41, 1241 (1961) [JETP 14, 886 (1962)]; Zh. 

Eksp. Teor. Fiz. 42, 1658 (1962) [JETP 14, 1151 (1962)]. 



CO cr 
— K- 



CO CT 
—* — 



CO CT 




CO -CT 
— » — 



CO CT 




CO -CT 
—4 



-CT 



n -cr 



CO CT 

4— 



CO CT 
-4 



(a) 



CT 



CO CT 
4— 



CO -CT 



(b) 



CT v w -CT 



co cr I -cr 

4 fr 



(c) 



Fig. 1 Diagrams contributes to the singularities 



9 



